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Abstract 

The relatively unexplored fast radio transient parameter space is known to be home to a va- 
riety of interesting sources, including pulsars, pulsar giant pulses and non-thermal emission from 
planetary magnetospheres. In addition, a variety of hypothesized but as-yet-unobserved phenom- 
ena, such as primordial black hole evaporation and prompt emission associated with coalescing 
massive objects have been suggested. The 2007 announcement by Lorimer et al. of the detection 
of a bright (30 Jy) radio pulse that was inferred to be of extragalactic origin and the subsequent 
consternation have demonstrated both the potential utility of bright radio pulses as probes of the 
interstellar medium and intergalactic medium, as well as the need for wide-field surveys character- 
izing the fast-transient parameter space. Here we present results from the 450 hour, 150 deg^ Fly's 
Eye survey for bright dispersed radio pulses at the Allen Telescope Array (ATA). The Fly's Eye 
spectrometer produces 128 channel power spectra over a 209 MHz bandwidth, centered at 1430 
MHz, on 44 independent signals paths originating with 30 independent ATA antennas. Data were 
dedispersed between and 2000 pc cm~'^ and searched for pulses with dispersion measures greater 
than 50 pc cm~^ between 625 fis and 5 s in duration. No pulses were detected in the survey, imply- 
ing a limiting rate of less than 2 sky~^ hour^^ for 10 millisecond duration pulses having apparent 
energy densities greater than 440 kJy ^s, or mean flux densities greater than 44 Jy. Here we present 
details of the instrument, experiment and observations, including a discussion of our results in light 
of other single pulse searches. 

1. Introduction 

The last decade has seen an explosion of interest in time domain radio astronomy. Driven 
by new wide field survey instruments, multi-beam receivers and computing advances, exploration 
of this regime presents opportunities to shed new light on known phenomena and perhaps reveal 
previously unseen processes as well 
can be conveniently divided into two categories: those that deal with slow transients, events lasting 
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from hours to years, and those that deal with fast transients, shorter duration events lasting from 
nanoseconds to seconds. Slow transients may originate from a diverse set of phenomena, ranging 
from the radio counterparts to relatively nearby supernovae or explosive events observed at very 
high redshift to variation in flux from radio jets associated with accretion disks, as discussed at 



length in Bower et al. (2007), Ofek et al. (2010), Bower & Saul (2011) and Becker et al. (2010). 



Fast radio transient phenomena have been observed to arise from only two fundamental physical 



origins: electrostatic discharge in planetary atmospheres, e.g. Zarka et al. (2008), and non-thermal 



magnetospheric emission from planets e.g. Zarka (2007) or stars e.g. Lorimer & Kramer (2005); 



Hallinan et al. (2008). By far the most commonly explored source of fast radio transient emis- 



sion are pulsars, including their relatively more intermittent flavors: "rotating radio transients" 
(RRATs), nulling pulsars and those producing giant pulses. A variety of as- yet-unobserved sources 
of fast transient emission have been proffered, among the more commonly cited are evaporating pri- 
mordial black holes ( Rees||1977 ) and compact object coalescence ( Hansen &: Lyutikov||2001 ). Other 
more speculative sources include emission from cosmic strings ( Vachaspati||2008 ) and beacons from 



extraterrestrial intelligence (Benford et al. 2008) 



The fast/slow division arises primarily from the different detection techniques employed in 
each case. Slow transients are identified at radio wavelengths using essentially the same techniques 
as used at optical wavelengths; images (usually synthesized from interferometric antenna arrays) 
from multiple epochs are differenced and thresholded, with the period between epochs driven by 
the parameters of a given experiment. Fast transients are generally identified via a method well 
established from decades of pulsar searches: high time resolution power spectra are corrected 
for the dispersive effects of the interstellar medium, dedispersed, assuming different integrated 
column densities of free electrons, dispersion measures or DMs, then collapsed to a time series and 
thresholded for impulsive events. In the case of periodicity searches a Fourier transform or fast 
folding algorithm is also applied to the time series. Although the existence of nulling pulsars has 



been known for some time (Backer 1970), only fairly recently have pulsar searches targeted them 



specifically. Following the McLaughlin et al. (2006) discovery of the extremely intermittent class of 



neutron stars dubbed RRATs, several archival data sets were reanalyzed including specific single 
pulse search algorithms in addition to periodicity searches. Although nearly all RRATs and nulling 
pulsars have been shown to possess an underlying normal periodicity, the duty cycle of emission can 
render some pulsars undetectable in a folded profile or Fourier transformed power spectra. For very 
low duty cycle objects or those with extreme variation in pulse flux, single pulse searches can be 
much more effective than periodicity searches and sometimes the only viable method of detection. 

In late 2007, a particularly perplexing radio transient was observed in an archival search of 
1.4 GHz pulsar survey data from the Parkes multibeam receiver ~ the so-called "Lorimer Burst" 



(Lorimer et al. 2007). This very bright single pulse was detected at a very high signal-to-noise 



ratio, with a peak flux nearly 100 times the search threshold, in a pointing a few degrees south 
of the Small Magellanic Cloud (SMC). The pulse clearly exhibited the quadratic chirp expected 
from an astrophysical pulse modifled by a cold plasma dispersion relation, with an inferred DM 
of 375 pc cm~^. Even liberal models for the galactic and SMC contribution to the total implied 
electron column density could account for only a fraction of the DM measured. Assuming the rest 
of the dispersion was due to a Milky Way-like host interstellar medium (ISM) contribution and 
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traversal of the much more rarified intergalactic medium (IGM), the lower Umit on the distance 
to the source was calculated to be ~600 Mpc. Suffice it to say, the implied energy release of 
~ 10^*^ ergs presented a challenge for astrophysical theory, and motivated wide ranging speculation 



of possible origins. Several subsequent searches (Keane et al. 2010 Deneva et al. 2009) did not 



detect any similar events, implying that such events must be exceedingly rare. Burke-Spolaor 



et al. (2010) presented the detection of several additional impulsive events in Parkes survey data 



with similar dispersive chirps to that seen in the Lorimer Burst but exhibiting clear indications of 



terrestrial origins. While the Burke-Spolaor et al. (2010) events showed an approximately quadratic 



frequency evolution, as would be expected for an astrophysical event, the detection of the events 
in multiple receiver beams simultaneously clearly points to a terrestrial source and the irregularity 
of received flux across the observing band and large pulse width differentiate them markedly from 
the Lorimer Burst. Recently, another possibly extragalactic burst was discovered in additional 



re-analysis of Parkes survey data (Keane et al. 2011), lending some support for the existence of a 



bonafide population of very bright extragalactic fast transient sources. Regardless of the source 
of such bursts, a population of extragalactic objects or events producing extraordinarily energetic 
radio pulses would provide an invaluable probe of the ionized IGM. 

Here we present a search using the 42-dish Allen Telescope Array for bright dispersed radio 
pulses, with specific attention paid to those of possible extragalactic origin. Section [2] describes 
the digital spectrometer developed for this experiment, installation, verification and calibration 
procedures. Sections [3] and |4] detail observations and analyses and Section [5] presents our results 
and interpretation. 



2. Instrument and Installation 

The single pulse search described here used the Allen Telescope Array (ATA) ( [Welch et al. 



2009) in an unconventional non-interferometric mode. Rather than pointing all 42 dishes in the 
same direction, each dish was pointed at a unique position, similar to a "fiy's eye." Such a mode 
yields a dramatically increased field of view at the expense of sensitivity, well matched to detecting 
bright, rare events. The primary half-power beam width (HPBW) of the ATA is approximately 2.5 
deg at 1.4 GHz, yielding a potential field of view of more than 200 deg^. 

In this experiment, each antenna signal path was processed independently using a purpose- 
built digital spectrometer. This device, dubbed the Fly's Eye Spectrometer, was constructed using 
the modular instrumentation infrastructure developed by the Center for Astronomy Signal Pro- 



cessing and Electronics Research (Werthimer et al. 2011). The full system consists of eleven field 



programmable gate-array (FPGA)-based 'iBOB' computing boards, each equipped with two 1024 
Msample/sec 'iADC' analog-to-digital converter cards. Each ADC board digitizes two independent 
single-polarization signal paths at 838.8608 Msamples/sec. The Nyquist sampled band is digitally 
down converted and decimated to a bandwidth of 209.7152 MHz within the FPGA and channelized 
using a 2^ channel complex biplex-pipelined polyphase filterbank. Power spectra are detected and 
accumulated for 625 fis, packetized into Ethernet UDP packets on each of the eleven FPGA boards 
and transmitted to a single Linux PC via an Ethernet switch. The payload of each UDP packet 
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contains a 21 byte header, which specifies a board ID, accumulation number and any error condi- 
tions, followed by 512 bytes of spectral data (unsigned byte power measurements x 128 frequency 



channels x 4 inputs). The entire Fly's Eye digital back end is described in detail in McMahon 



(2008) and Siemion et al. (2010). Large portions of the hardware design and software are open 



source and freely available at http://casper.berkeley.edu. 

During the course of analyzing our initial observations, a subtle timing error was discovered in 
the accumulation counter produced by the Fly's Eye hardware. This error was initially attributed 
to a flaw in the digital design for the instrument involving the 1 pulse per second synchronization 
logic. As a result, we disconnected the 1 pulse-per-second input signal from each iBOB board 
prior to beginning our observation campaign. It was later discovered that the principal cause of 
the timing error was not the flaw in the digital design, but rather that the sampling clock was 
not properly locked to the observatory reference. By comparing the unix time stamps applied to 
each spectra by our data collection computer with the hardware counter applied by the Fly's Eye 
spectrometer, we determined that the sampling clock frequency was only absolute to about 1 part 
m lO'^. This level of error translates to a center RF frequency and bandwidth ambiguity of about 
84 kHz. Because the Fly's Eye spectrometer derives integration time from counting the sampling 
clock, the unlocked clock synthesizer also imposes an ambiguity in integration time of about 60 
nanoseconds. In total, these effects correspond to an additional ~0.03% temporal smearing at 
1430 MHz. Although a small amount of additional temporal smearing has only a small impact on 
single pulse searches, it renders simple barycentering and folding of known pulsars impossible. We 
accounted for this effect during operability determination by 'searching' for a known pulsar rather 
than directly folding on its known ephemeris. 

The Fly's Eye Spectrometer system was installed at the Allen Telescope Array in December 
of 2007 and connected to the 44 antenna-polarizations exhibiting the lowest system temperature. 
Galactic HI detection along several lines of sight served as a coarse operability check. To further 
test the Fly's Eye Spectrometer, and most importantly to test our ability to detect bright dispersed 
pulses, we observed several bright pulsars including PSR B0329+54, the brightest known pulsar in 
the northern celestial sphere in terms of mean flux and also PSR B0531+21 (the Crab pulsar) - a 
canonical source of bright dispersed radio pulses (Bhat et al. 2008 Cordes et al. 2004). Figure [l] 



shows folded pulse profiles for a 1290s B0329+54 integration as observed in each of the 44 Fly's 
Eye inputs. All inputs but one (8B) show some detection, with only 4 others (9D, 6C, 4D, 2C) 
showing significantly degraded signal-to-noise. Figure [2] shows an unweighted incoherent sum of 
the same data plotted alongside a reference profile. 



The signal-to-noise ratio of a detected pulsar pulse profile can be calculated as (Cordes & 



Chernoff|[l998| ): 

^/N=^Vn,t,n,AfN^ (1) 

»-'sys 

where ^/jy is the detected signal to noise ratio, ^niean is the mean flux density of the pulsar, 5sys is 
the system equivalent flux density (SEED) of the observing system, Up the number of polarizations 
summed, tint the integration time. A/ the bandwidth observed and N^^ is the number of harmonics 
used in a harmonic sum, which depends on the pulse period P and pulse width W as A^^^^ ~ P/W. 
For the case of multiple antennas sampled synchronously and detected independently, rip is equal 



- 5 - 



S 6 




o 

°3 5 





ATA Fly's Eye 
PSR B0329+54 - 12/2007 






























V 






B C 

Signal Path 



Fig. 1. — : Folded pulse profiles of PSR B0329+54 as observed in individual Fly's Eye Spectrometer 
inputs for a 1290s integration at 1430 MHz. Two full turns plotted for clarity 



to the total number of antenna-polarization signals incoherently summed, nant-pol- Applying this 
equation to the summed profile shown in Figure [2] with n-ant-poi = 44 and the expected mean 
SEED for individual ATA antennas (~10 kJy), we infer a flux density for B0329+54 of Smean ~ 100 



mJy. This value is in reasonable agreement with a mean value extrapolated from Manchester et al 
(2005), Suiean = 190 mJy, especially considering that the pulsar B0329+54 has a variable mean 



observed flux density of a factor of ~3 at 1400 MHz (Wang et al. 2005: Liu et al. 2006). PSR 



B0329+54 is thus a poor flux calibrator, but it is never-the-less a useful diagnostic source for a low 
sensitivity pulse search. 

Observation of giant pulses (GPs) from the Crab pulsar provided the final end-to-end test of 
the Fly's Eye observing system. Crab GPs are known to follow a power law brightness distribution, 
often parameterized in terms of a cumulative probability distribution P{Ei > Eq) = KEq, where P 
gives the probability of a pulse having a pulse area Ei greater than Eq. Here pulse area is defined as 
Ei = SiWi, where Si is the mean intrinsic flux density of the pulse over an intrinsic time Wi. While 
other authors have referred to the quantity E[ as "energy" , in later portions of this work we will 
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Fig. 2. — : PSR B0329+54 folded pulse profile as detected by the incoherent sum of all 44 inputs to 
the Fly's Eye Spectrometer (shown in[T]) for a 1290s integration at 1430 MHz (left) and a reference 



profile from Gould & Lyne (1998) taken at 1408 MHz 



use the slightly more accurate "energy density." At 1300 MHz, Bhat et al. (2008) gives a ~ -1.9 
for energy densities greater than 10 kJy /xs, with K = 4.7 x 10~^ and Eq in kJy //s. Rearranging 



equation 3 in Deneva et al. (2009) gives an expression for the minimum detectable energy density, 



mS< 



sys 



W 



(2) 



W is the observed pulse width, usually taken to be the quadrature sum of the various sources of 
broadening, both astrophysical and instrumental and m the signal-to-noise threshold. For a pulse 
with an observed width limited by our digital hardware and assuming a single polarization SEFD 
of ~10 kJy, a characteristic m = ba minimum detectable energy density for the incoherent sum of 
19 inputs is 20 kJy /is or an mean flux density of ~32 Jy for a 625 /US pulse. Based on the expected 
Crab GP distribution, we should observe a pulse with Ei > 20 kJy ^us about 18 times per hour 
on average. Figure [s] shows the detection of 10 bright GPs at the expected DM of 56.8 pc cm~^ 
for the Crab pulsar in a 50 minute unweighted incoherently summed observation using the 19 best 
performing FE inputs, as determined from Figure [T| 

The set of 44 antenna inputs ultimately used for the FE observing campaign originated with 30 
independent antennas, 14 of which included both X and Y polarizations. Figure |4] shows antenna 
performance, described by the system equivalent flux-density (SEFD), for each of the 44 inputs 
chosen for the Fly's Eye observing campaign. The values given here were determined by interfer- 
ometric observation of standard calibrators interleaved with Fly's Eye observation. Details of the 
Fly's Eye instrument parameters are given in Table [l| 



3. Observations 

During the period February 2008 to May 2008, we conducted approximately 480 hours of 
drift-scan observations with the Fly's Eye Spectrometer ( Table [2]). Data were collected in 60 
minute intervals, each consisting of 58 minutes of drift observation followed by 2 minute diagnostic 
observations used for monitoring the health of the telescope and instrumentation. Fly's Eye obser- 
vations produced data at a rate of roughly 36 GB / hour, resulting in approximately 18 TB total 
collected data for the entire observation period. The data are archived in Berkeley and available 
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Instrument: Fly's Eye 



MJD: 54456.418 
F : 1431.6 MHz 




10 15 
Sigrml-Lo 



EL 

20 




Max SNR: 23.8 
E Inputs: 19 



50 100 
DM (pc or 



1500 2000 
Time (s) 



Fig. 3. — : A standard single pulse detection plot for a ~50 minute observation of the Crab pulsar 
after summing 19 Fly's Eye inputs. Several giant pulses are apparent at a DM of 56.8 pc cm~^. 
Other features are radio frequency interference. From top, left to right, the panels show a histogram 
of detection signal-to-noise ratio (SNR) for pulses > 5a, a histogram of detection dispersion measure, 
detection dispersion measure against detection SNR and detection time vs. dispersion measure 
with detection signal to noise indicated by plot point radius. This plot produced using PRESTO 



(Ransom 2001) 



for analysis by request to the authors. At the time of these observations, the ATA was undergoing 
commissioning, and a variety of system performance issues were being actively addressed. The 
variation in SEED from antenna to antenna and less-than-complete utilization of the 42 installed 
antennas are reflective of these issues. To aid in dynamically determining signal path operability, 
a fixed pointing strip along a constant declination angle of +54°, in which antennae were spaced 
1 half-power beam width apart, was chosen for drift scan observations. As the bright pulsar PSR 
B0329-I-54 drifted through the beam pattern at the sidereal rate, its detection or non-detection 
was used to determine whether or not a particular signal path was operable. Declination +54° is 
well away from any significant source of galactic electron density confusion, the median maximum 
galactic DM contribution along this path is ~68 pc cm~^ (from the NE2001 model, Cordes & 



Lazio (2002)). Figure [S] illustrates the overall observing efficiency after applying the B0329+54 
detectability metric. Out of a total of 921.1 input-days of observing, 579.9 input-days showed the 
expected detections of B0329+54, for a total observing efficiency of ~63%. The poor efficiency 
in Epoch 1 likely reflects an unaccounted for change in our pointing script that directed antennas 
away from +54°. Although we believe most signal paths were operable, we have conservatively 
excluded these data. 
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Fig. 4. — : Mean antenna-polarization performance for each input to the Fly's Eye Spectrometer, 
as determined by interferometric observation of standard calibrators interleaved with Fly's Eye 
observations between 02/2008 and 05/2008. Errors are it la. 



4. Analysis 
4.1. Data Preparation 

Power spectra time series for each of the 44 inputs to the Fly's Eye Spectrometer were extracted 



as individual "filterbank" format files (Lorimer et al. 2000), broken into analysis chunks of length 
2^*^ samples (representing 12 minutes). This length was chosen to allow an entire analysis chunk and 
set of dedispersed time series to be kept in computer memory during analysis. Prior to dedispersion, 
power spectra were normalized or "equalized" across both frequency channel and time. Equalizing 
across frequency channels has the primary effect of correcting for the rippled bandpass imposed 
by both analog filter response and digital down conversion. Equalization of accumulation values 
mitigates broadband gain changes and broadband impulsive interference. This process was carried 
out as follows. 



We denote the power in channel i at (discrete) time t as Pi{t) £ [0,255], and we computed a 
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Table 1:: Fly's Eye Instrument Parameters 



Center Frequency Vo 1430 MHz 

Number of Channels A^chan 128 

Channel Width Avi 1.6 MHz 

Bandwidth 210 MHz 



Beam Width 9 2.5 de^ 

(per antenna, HPBW) 

Solid Angle Q. 147.3 deg^ 
(Instantaneous) 

System Temperatur^ T^ys 50 K 

Gain'^ G 6.25 x 10"^ K/Jy 

Dish Diameter D ~ 6 m 



MacMahon & Wright ( 2009 1 



''Assuming all signal paths are operable but 
accounting for some dual-polarization olaservations, see 
Section |2] 

"^Nominal value 

Table 2:: Fly's Eye Observations 02/2008-05/2008 



Epoch MJD Efficiency 

1 54512.24 - 54515.97 5% 

2 54519.20 - 54521.68 50% 

3 54526.13 - 54528.64 86% 

4 54533.04 - 54535.63 73% 

5 54540.17 - 54542.59 75% 

6 54547.18 - 54549.60 86% 

7 54554.18 - 54556.62 93% 

8 54561.26 - 54563.62 66% 



mean power per channel 



^0 t=0 



over some time period Tq. Tq is typically set to the length of an analysis chunk, 2^*^ samples. A 
particular value Pi{t) is then divided by the mean Pi. i.e. the equalized value -P/(t) = Pi{t)/Pi. 
The mean power in each channel is then unity, since 



^^£p/(t) = l (4) 



t=0 
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Fig. 5. — : Diagram showing operability of each antenna/input as a function of epoch, based on 
detectabihty of B0329+54. Filled circles indicate that a given signal path is operable. For an 
input/epoch pair to be considered operable, B0329+54 must have been detected at every opportu- 
nity within the epoch. On the Y axis are each spectrometer input, labeled by their ATA antenna 
identifier followed by spectrometer input number. The suffix on the antenna identifier indicates 
which of two dual linear polarization feeds was used. 

Mean power equalization was performed on the frequency spectrum equalized values -P/(t). The 
power mean over all frequency channels for a single integration (time sample t) is defined as 

i=0 

N is the number of channels. With the mean powers P'{t), we can define the equalization of the 
powers P'i{t). The mean power equalized values P'i'it) = P[{t)/P'{t). This procedure ensures that 
the sum of the power samples for any time Tq is normalized such that X]£o^ ^ii^) — effectively 
flattening the DM = time series. 

Prior to the equalization process, individual frequency channels with especially large amounts 
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of interference were identified and logged. Our algorithm used the variance of each frequency 
channel over an analysis chunk as a measure of the amount of interference in that channel. Using 
the previously defined quantities, we computed the variance of the values Pi{t) for each channel i 
over To, and then fit a polynomial to the resulting curve Var(Pj) (Figure [g]). Frequency channels 
for which the computed variance differed from the polynomial fit by Var(Pj) > 2a were excluded 
from subsequent dedispersion. We explicitly excluded 8 frequency channels at the top of the band 
and 13 frequency channels at the bottom of the band due to analog filter roll off and the presence 
of bright air route surveillance radar below 1350 MHz. 




3 2 64 96 128 



Frequency Channel 

Fig. 6. — : Variance vs. frequency channel for a 2^^ spectra Fly's Eye observation (red), an iterative 
polynomial fit to the variance curve (green) and ±30" bounds on the fit (black). Band edges are 
excluded a priori from the polynomial fit, as instrumental response is poor in these regions due to 
filter roll-off. Band edges and individual frequency channels with high variance are excluded from 
the de-dispersion process. 



4.2. Single Pulse Search 



A single pulse search in dedispersed time series for events having a signal-to-noise ratio (SNR) 
greater than five standard deviations above the mean, a > 5.0, was carried out over 744 trial 
DMs between 0-2002 pc cm~'^ using the SigProc tools (Lorimer et al. 2000). An approximate 



matched filtering algorithm was employed to increase sensitivity to broadened pulses in which each 
dedispersed time series was iteratively smoothed by adding 2" adjacent time samples over the range 
n = to 10, following Cordes & McLaughlin (2003). This results in an effective box car smoothing 



of maximal window size 2^*^ samples or 0.64 to 5.12 s, depending on the level of time collapse (Table 
[3|. The initial DM step size was set such that the time delay associated with the DM step was 
less than the sampling interval. Spectra were iteratively collapsed in time by a factor of two at 
DMs 329.0, 658.0, 1314.0 pc cm^^ to speed analysis, as detailed Table [sj At these thresholds a 
pulse would subtend a minimum of 2, 4 or 8 spectra at the top of the band, respectively, and thus 
halving the effective time resolution imposes no loss in sensitivity. Similarly, the DM step size can 
be reduced by the same factor and remain less than the new effective time resolution. These search 
parameter choices were all designed to ensure that the dominant source of temporal smearing was 
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due to the unavoidable (at the time of analysis) integration time smearing and in-band dispersive 
smearing. Like many problems in radio astronomy, the analysis of multibeam pulse search data is 
readily parallelizable. Here we distributed 58 minute observations to individual compute nodes, and 
parallelized inputs over individual CPU cores. In total, our analysis consumed ~20,000 CPU-hours. 

Table 3:: Single Pulse Search Parameters 



Dispersion Measure Range ADM Time Collapse Factor 

0.0 - 329.0 pc cm"^ 1.0 1 

330.0 - 656.0 2.0 2 

658.0 - 1310.0 4.0 4 

1314.0 - 2002.0 8.0 8 



4.3. Post Processing RFI Filter 

After dedispersion, any strong signal present in a dynamic spectra - be it from interference 
or a real event - will be detected at multiple DMs, strengths and times. The distribution of 
these detections depends on the observed properties of the signal, the range of DMs searched and 
any pre- or post-dispersion processing applied. In the case of quadratically chirped radio pulses, 



this distribution follows a characteristic functional form (Cordes & McLaughlin 2003). Likewise, 
certain kinds of interference will exhibit predictable detection distributions. Wideband, temporally 
narrow RFI will be detected over many DMs with the highest strength detections at low DMs. 
Narrowband, long duration RFI will also be detected at many DMs but will peak in strength when 
the dispersion path for a trial dispersion measure optimally overlaps the narrow-band interference. 
Finally, wideband and long duration interference or rapid gain changes will cause an excess of 
detections in all DMs for the duration of the event. 

As a first cut on the vast number of high SNR candidates detected, we flagged any 1 second 
input-time region where the highest SNR candidate in that region was < 50 pc cm~^ or > 1950 pc 
cm~^, or the total number of pulses over 5cr in the same DM regimes exceeded a factor of 4 times 
the mean number of pulses in each regime. The DM < 50 pc cm~^ would reject any relatively 
nearby galactic events, but distinguishing true astrophysical bursts from interference at these low 
DMs is very difficult because of the correspondingly small quadratic chirp. Further, our focus 
was primarily on potentially extragalactic events for which the DM contribution from the Milky 
Way and host ISM should well exceed this threshold. Figure [7| shows the results of applying this 
metric to pulse detections from two observations of the Crab Pulsar. The algorithm was effective 
at rejecting strong interference and avoided rejecting true astrophysical events. Strong interference 



dominates Figure 7a , but is greatly diminished in 7a Figures 7b and |7d| shows the detection of a 



bright pulse left untouched after applying our algorithm. 

Following RFI rejection, all events with a SNR o > 8.0 detected in 2^*^ sample analysis chunks 
for which the mean SNR of all detections was a < 5.5 were selected for visual analysis. Plots similar 
to Figure [7| were examined for each detection. If an event did not appear to be interference, we 
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Crab Pulsar Observation: Single Pulse Detections (Raw) 
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Fig. 7. — : Detection of giant pulses from the Crab pulsar for two individual FE inputs before (7a 



7b) and after (7c 7d) application of a post processing RFI filter. Each plot shows events vs. time 
and trial dispersion measure, here the radii of plot points are proportional to (signal to noise)^ 
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extracted and closely examined t vs. u spectrograms of 1, 2, 4 and 10 seconds around the event. 
Figure [8] shows the SNR distributions for all pulses detected in operable signal paths, all pulses in 
operable paths after applying the first cut RFI rejection algorithm described above, and (inset) all 
pulses detected in observations having a mean pulse detection SNR a < 5.5. Upon close inspection, 
none of the pulse candidates identified appeared to be of astronomical origin. An example of the 
pathological interference that escaped our interference rejection algorithms is shown in Figure |9j 
Strong pulses were detected in regions where dedispersion curves aligned with the triangle-wave 
modulation of the interferer, at a DM of ~80 pc cm^'^. This particular interferer was detected at 
multiple epochs, and appears to originate with an orbiting satellite. 
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Fig. 8. — : SNR histograms for all pulses detected in operable signal paths, all pulses in operable 
paths after applying the first cut RFI rejection, and (inset) all pulses detected in 2^*^ sample 
observation chunks having a mean pulse detection SNR a < 5.5. 



5. Discussion 

Our results indicate that the millisecond radio sky is relatively quiescent at the energies probed 
by our experiment. For a threshold of 8a, our minimum detectable energy density is -Emin = 111 
kJy /is or a mean flux density of ~178 Jy for a 625 fis pulse. Based on our non-detection, pulses 
of these energy densities originating with an isotropic progenitor must occur at a rate less than 2 
sky~^ day~^. We can determine the limiting rate of occurrence of bursts as a function of antenna 
sensitivity by computing the rate 

T^iE > E^^) = j2n.TMSEFD>SEFD^,^) 
where -Bmin corresponds to the detectable limit for pulses of a given duration in an antenna having 



SEFD^iji. Figure 11 shows the event rate limit calculated from these results alongside other 
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Fig. 9. — : A spectrogram plot of intensity in the time-frequency plane of a satellite interferer that 
passed our RFI rejection algorithms. Edges of the observed band, where system response is poor, 
have been masked. 



recent single pulse searches in Table |4j Our rate limit does not yet sample potential coherent radio 
emission processes from neutron star binary inspirals or gamma ray bursts, but we can place an 
order of magnitude limit on the luminosity of coherent radio emission from core collapse supernovae 
(CC SNe). Assuming isotropic emission, our survey could have detected ~20 10 ms events with 
an intrinsic rate of 1000 sky~^ day~^ at an apparent energy density limit of -Eumit = 10^ kJy fis. 
Assuming that coherent radio emission from CC SNe is beamed over a solid angle 17, we can limit 
the emission cone to 17 < 47r/20. Translating this into an upper limit on luminosity using the radius 
of a spherical volume of 1 Gpc'^ , we have 



47ri?2j^^ 1 



20 At 

For emission over a bandwidth Av = 1 GHz and At = 10 ms we have L < ~ 2 x 10^^ erg sec~^. 



(7) 



As discussed in Lorimer et al. (2007) and Keane et al. (2011 ), if the two detected pulses inferred 



to be at cosmological distances are indeed real, they must represent an entirely new source class. 
The extreme SNR of the two sparse detections is curious, and as previously discussed by several 
authors, contradicts the assumption of an isotropically distributed population. In the case of this 
experiment, we would have expected to detect ~15 events similar to the [Lorimer et ah (2007) 
burst, but <C 1 similar to the Keane et al. (2011) extragalactic event. The incompatibility of 
the implied rates for the two isolated detections is difficult to explain. The galactic latitude of the 
two events differ significantly, b = —41.8° for the Lorimer et al. (2007) event and b = —4° for the 



Keane et al. (2011) burst. While one might guess that being much closer to the galactic plane 



would make it easier to explain a large dispersion measure, of the 227 known pulsars between 3° 
and 5° off the galactic plane, none has a DM > 460 pc cm~^, standing in stark contrast with the 



^Assuming an Euclidean isotropic distribution and correcting for the Fly's Eye survey's lower sensitivity. See, for 
example, Deneva et al. (20091 for an exposition of this calculation. 
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Keane et al. (2011 ) event at a DM of 745 pc cm ^ However, the increased intragalactic path length 



of this detection does present more opportunity for an miseen highly ionized nebula to make an 
aberrant contribution to the total integrated free electron column density. 



The 15 Burke-Spolaor et al. (2010) detections also remain puzzling. Assuming the phenomena 



that generated these bursts is not unique to the Parkes site, we can estimate the number of similar 
events that the Fly's Eye survey should have detected. Assuming that these bright events were 
observed in far out side lobes, a reasonable comparison between the two surveys reduces to simply 
the ratio of their observing time, thus the Fly's Eye survey should have detected 



450 , 

15 events x h ~ 10 events (8) 

346.1 + 532.4 ^ ^ 

of these events as well. Note that here we assume antenna efficiency and system temperature 

differences are negligible. We will refrain from speculating on the cause of these events, but note 

that our observations well sampled the diurnal cycle as well as varying levels of precipitation. We 

are unaware of any lightning activity in the near vicinity of the ATA during our observations. 

It remains perplexing as to why all of these unique transient bursts have been detected only at 
the Parkes telescope. However, with the aggressiveness with which interference must be excised in 
such experiments, and the varying means by which it is accomplished, we speculate that it is not 
wholly out of the question that some processing pipelines could be better tuned to detecting single 
pulses at extragalactic DMs or other unexpected characteristics. The fact that all of these pulses 
were themselves discovered in reanalyses of previously mined surveys, with the second extragalactic 
event discovered in a 3rd reanalysis, indicates that other extant surveys may harbor additional 
as- yet undetected events. Our own pipeline has been honed by the exercise of this experiment, 
and a future search would undoubtedly be superior. For instance, summing polarizations would 
yield an additional \/2 sensitivity on some antennas. We attempted multi-beam coincidence RFI 
rejection, excluding pulses found at similar DMs and similar times in multiple antennas, and found 
it somewhat ineffective. We attribute the poor performance of this algorithm to a number of factors, 
among them that our algorithm didn't account for variations in signal path sensitivity and multi- 
path propagation of interferers led to strong variation in detection signal-to-noise from antenna to 
antenna. 

In previous searches where multi-beam excision has been effective, multiple receivers were 
co-located in the focal plane of a single telescope. Our results indicate that the interference en- 
vironment is generally better correlated in these more confined configurations. Even so, applying 
a multi-beam coincidence metric in the time/frequency plane, prior to baseline subtraction and 
thresholding, would undoubtedly perform better. We explored these possibilities in parallel with 



the analysis described here, in anticipation of a future reanalysis, in Hogden et al. (2011). Of the 



methods explored on a subset of Fly's Eye data in Hogden et al. (2011), we find that a combina- 



tion of Huber filtering and adaptive interference cancellation performs optimally for the types of 



interference seen in this experiment. Friends-of- friends algorithms, as described in Deneva et al. 
(2009) or Hough transform-based detection ( Fridman|2010 ) would also improve our ability to detect 
low SNR events. In any case, additional surveys and analyses will soon detect more examples of 
isolated radio pulses which will expand the available sample. New wide-field interferometric radio 
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telescopes will be well equipped to not only detect additional events, but also provide much better 
localization than has so far been possible with single dish facilities. Precise localization would be 
extremely useful in differentiating between true extragalactic events and those that may appear so 
because of passage through highly ionized intragalactic regions. 



6. Summary 



We have developed a novel multiple-input digital spectrometer which we have used to conduct 
a wide field search for bright dispersed radio pulses using the Allen Telescope Array. This wide 
field search yielded no detections, allowing us to place a limiting rate of less than 2 sky~^ hour~^ 
for 10 millisecond duration pulses having mean apparent flux densities greater than 44 Jy. The 
flux densities probed by this experiment are well above individual pulses from known pulsars and 
RRATs, just grazing the very brightest of the giant pulse producing pulsars, none of which are 
present in the field surveyed. We have placed new limits on very bright coherent emission from 



events similar to the singular event described in Lorimer et al. (2007). Our results indicate that 



the Lorimer et al. (2007) event must belong to a very rare source class, if it is indeed astrophysical. 



We did not detect any quadratically dispersed terrestrial interference similar to that seen at the 



Parkes observatory, e.g. Burke-Spolaor et al. (2010), consistent with other non-Parkes surveys, e.g. 



Deneva et al. (2009). 



The work presented here has shown that sources of bright fast transient radio emission must 
be relatively rare. We have also shown both the utility, and associated of challenges, of using an 
interferometric array in a multi-pointing "fly's eye" mode. The next generation of radio inter- 



ferometers currently being built in preparation for the Square Kilometer Array (Dewdney et al 



2009 ) will offer new opportunities to explore the fast transient regime with greater sensitivity over 



large solid angles. The use of these new instruments for fly's eye mode surveys will require mak- 
ing individual antenna or station data accessible to sufficient digital hardware, and we encourage 
this consideration to be taken into account during the design phase. Commensal surveys for fast 
transient emission with interferometers, using the incoherent sum of antennas pointed in the same 
direction, could offer an excellent trade-off between sensitivity and solid angle while incurring little 



additional hardware cost and no additional observing time, see e.g. Macquart (2011). Again, such 



capabilities will require consideration early in the design process in order to be realized efficiently. 
The exploration of the fast radio transient parameter space is just beginning, and the contrasting 
results of this and other experiments clearly indicate we have much yet to learn. 
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Fig. 10. — : A plot of system sensitivity, described by the system equivalent flux density, vs. total 
operable observing, described by the observing solid angle • time product. Plot includes both 
operable observing, taken from the set of signal paths where PSR B0329+54 was consistently 
detected, and the final low-RFI analysis set. In cases where two polarizations were observed for a 
given antenna, the lower SEFD was used. 
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Fig. 11. — : Pulse energy density vs. rate limit for the surveys in Table |4j Rate limit curves assume 
a 10 millisecond pulse duration. Shaded bars on the Lorimer Burst and the extragalactic event 



described in Keane et al. (2011) represent 2 sigma confidence (Gehrels 1986) on the Euclidean 



isotropic distribution. The point for terrestrial events identified in Burke-Spolaor et al. (2010) 



assumes off-axis detection. Rates of core collapse supernovae, gamma ray bursts and binary neutron 
star inspirals in a 1 Gpc^ volume are taken from Madau et al. (1998), Guetta della Valle (2007) 



and Kalogera et al. (2004), respectively, via Lorimer et al. (2007). Here we assume flat sensitivity 



across the HPBW of each receiver beam and do not take into account reduced sensitivity in wide 



field-of-view sidelobes, except in the case of this work and Deneva et al. (2009). The curved line 



for the Siemion et al survey reflects variation in antenna system temperature. 
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Table 4:: Recent L-band Single Pulse Surveys 



Survey 



Tobs Total Solid Angle E„iii|^ 

(hours) (dog'^) kjy /is 



Edwards et al. ( 2001 1 



Re-analysis by 



Burke-Spolaor & Bailes 



(20101 



Manchester et al. 



Re-analysis by 



(20061, 



Lorimer et al. 



(20071 



Jacoby et al. 



(20091 



Re-analysis by Burke-Spolaor & Bailes 



(20101 



Manchester et al. 



Re-analysis by 



(20011, 



Keane et al. 



(20101 and 



Keane et al. 



(20111 



Deneva et al. 



(20091 



Fly's Eye (this work) 



346.1 

480.7 

532.4 

1864.3 
461.0 
136 



0.55 

0.556 

0.556 

0.556 
0.0187[^ 



147 



0.8 

0.8 

0.8 
0.1 

440 



"For a 10 millisecond pulse. Values given here are conservative, see references for details. 
'Parkes Multibeam - 14' HPBW/beam x 13 beams 
'^Arecibo ALFA - 3.5' HPBW/beam x 7 beams 
■^2.5° HPBW/beam x 30 beams 
•^Using a nominal SEFD of 8 kJy 
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